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Gluon multiplicity in coherent diffraction of onium on a heavy nucleus. 
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We derive the cross section for the diffractive gluon production in high energy onium- 
nucleus collisions that includes the low- a; evolution effects in the rapidity interval between 
the onium and the produced gluon and in the rapidity interval between the gluon and the 
target nucleus. We analyze our result in two limiting cases: when the onium size is much 
smaller than the saturation scale and when its size is much larger than the saturation scale. 
In the later case the gluon multiplicity is very small in the quasi-classical case and increases 
when the low- a; evolution effects in onium become significant. We discuss the implications 
of our result for the RHIC, LHC and EIC phenomenology. 

I. INTRODUCTION 



Diffractive gluon production in high energy pA collisions and in Deep Inelastic Scattering (DIS) 
is a sensitive probe of the Color Glass Condensate characterized by 

high parton density and gluon saturation [ill , [ijj]. Diffractive processes played a pivotal role in 



identif ying the first signatures of the gluon saturation in DIS at HERA 



19, 3 



P, Q, Q, 



16, 



18, 



2l| . They are of great interest as a tool for studying the low-x dynamics in pA collisions 



at RHIC and LHC as well as in DIS at the proposed EIC. Study of high parton densities in 
deuteron - gold collisions at RHIC has provided many novel insights into the structure of nuclear 
matter and has been focused on inclusive processes. By triggering on hadrons in the deuteron 
fragmentation region ("forward" rapidity) one is able to access very low values of Bjorken x that 
are sensitive to the gluon saturation. Investigation of energy, rapidity, centrality, and transverse 
momentum dependence of various production channels offers an opportunity to attain a better 
understanding of the nuclear and hadron structure at low x. Among the channels which have been 



discussed in t 
gluons 



quarks 



26, 
43, 



lis context are total hadron multiplicities 



27 



22 



28 



29 



30, 



31 



32, 



441 ] , prompt photons [45[] , di-leptons 



33 



34 



reavy 



46. 



47 



quarks 



23 
36 



37, 



38 



incusive production of 



39, 



40, 



42J, valence 



n 

481 ] and identified hadrons [49] (for a review 



2 



see e.g. |50l . |51||). Diffractive production in pA collisions offers another avenue for exploring the 
low-x dynamics. Motivated by a possibility to measure the diffractive production in pA collisions 
at RHIC and LHC we analyze in this paper diffractive gluon production in onium-heavy nucleus 
collisions. We intentionally avoid discussing the "dipole content" of the proton light-cone "wave 
function" and concentrate entirely on quantities that can be calculated in perturbation theory. Our 
results can be equally well applied to diffractive gluon production in DIS in which the light-cone 
"wave function" of the virtual photon is well-known. Diffractive gluo n production in DIS has been 



discussed in many publications 







56 



53, 



58 



59j ] and has been limited to the 



quasi-classical approximation and/or phenomenological models. In this paper we go beyond the 
quasi-classical approximation and include the low-x evolution effects at all rapidity intervals. 

The paper is structured as follows. In Sec. HI 
in the quasi-classical approximation derived in 



we review the result for diffractive gluon production 
55, 



60, 



611 ] . In Sec. IIIII we generalize these results by 



including the effect of quantum evolution. We consider separately the case when the rapidity gap Yq 
between the produced gluon and the target equals the gluon's rapidity y (Fig. [2]) and a more general 
case when y > Yq (Fig. [3|). The corresponding cross sections are given by (]14p and (|16p in terms 
of the dipole distribution in proton n p (r, r ; , b, y), the forward dipole-nucleus scattering amplitude 
N(r,b,y) and diffractive dipole-nucleus scattering amplitude N^(r, b, y; Yq). In Sec. HVI we review 
the main properties of n p (r, r', b, y) and iV(r, b, y) in the linear regime and in the saturation regime 
and demarcate the kinematic landscape. We then proceed in Sec. [V]by performing analysis of the 
diffractive gluon production in the quasi-classical approximation in various kinematic regions. The 
results are displayed in (|5ip and (|56f) . In Sec. IVII we do similar analysis in the case of low-x 
evolution, see (162p and (|65p . We summarize and discuss the phenomenological importance of the 
obtained results in Sec. IVII1 



II. DIFFRACTIVE GLUON PRODUCTION IN THE QUASI-CLASSICAL 

APPROXIMATION 

The quasi-classical approximation of the hadron-nucleus interactions is valid when a typical 
parton coherence length l c is much larger than the nuclear size Ra in the nucleus rest frame. 
The former is approximately given by l c ~ 2m N x , where m n is a nucleon mass. It follows that 
the quasi-classical approximation holds for x < R,Am N • Owing to the large coherence length, 
the process of diffractive production can be considered as proceeding in two stages: gluon emission 
long time before the collision followed by the instantaneous interaction. This picture is particularly 



FIG. 1: One of the diagrams contributing to the diffractive gluon production in onium (P) - heavy nucleus 
(A) collisions in the quasi-classical approximation. Notations are explained in the text. 



simple in the transverse configuration space since the parton transverse coordinates do not change 
in the course of instantaneous interaction. As the result, the cross section can be represented 
as a convolution of the proton's light-cone "wave-function" and the scattering amplitude in the 
transverse configuration space, see Fig.[TJ In the quasi-classical approximation, the cross section for 
the diffractive gluon production in onium-heavy nucleus collisions has been derived in 



55 



60 



61|. 



Using notations of 

da(k,y) ol s Cf 1 



551 ] . see Fig. [U it reads: 

2;, ,,2. ,2. I Z l ~ X 



d 2 k dy 



vr 2 (2vr) 2 



d bd z\d Z2 



X e - ik ( Z l- Z 2) / e --P( X :y: Z l) 



Zl — x| 

(x-y) 2 



Z2 



z 2 



|Z 2 



e -P(x,y,z 2 ) _ e -TW~ c 



|Z 2 

(x-y) 2 



(1) 



where x and y are the transverse coordinates of quark and antiquark, zi, z 2 are the transverse 
coordinates of the gluon in the amplitude and the complex-conjugate amplitude correspondingly, 
see Fig. [TJ The qqg propagator reads 



37 



61 



62] 



<-*P{ - p ! x - Y- z) I = ^T> ( -g(x - z) 2 Q 2 s0 - i(y - z) 2 Q 2 s0 + -L ( x - y) 2 Q 2 



The gluon saturation scale is given by 



4TT 2 a s N r . 



P T(h)xG(x,l/r 2 ), 



(2) 



(3) 



Ni-1 

where p is the nuclear density, T(b) is the nuclear thickness function as a function of the impact 
parameter b. The gluon distribution function reads 



xG(x,l/r 2 ) 



0! s Cf 



In 



1 



(4) 



vr r 2 A 2 ' 

with A being some non-perturbative momentum scale characterizing the nucleon's wave function. 
In the framework of the dipole model 63f] , the gluon evolution is easily accounted for in the 



large N c approximation. Let us introduce the forward elastic dipole-nucleus scattering amplitude 



4 



N(r, b, Y). In the quasi-classical approximation it reads 64| 



.1-2 r>2 



iV(r,b,0) = 1 -e"8 r . (5) 



At large N c , (HJ, © and © yield 



d^kdy n z (27r) 7 

zi - x zi - y \ / z 2 - x z 2 - y 



|zi-x| 2 |zi-y| 2 y \Jz 2 -x| 2 |z 2 -y| 2 
x [iV(x - y, b, 0) - AT(x - z 1; b, 0) - A(y - z ls b, 0) + A(x - z u b, 0)A(y - z u b, 0)] 

x [JV(x - y, b, 0) - 7V(x - z 2 , b, 0) - N(y - z 2 , b, 0) + JV(x - z 2 , b, 0)JV(y - z 2 , b, 0)] . (6) 

Integrating (|6j) over all transverse momenta yields delta function (2-7r) 2 5(zi — z 2 ). Hence, the 
total cross section per unit rapidity reads after a simple calculation 

da(y)_a s C F f (x - y) 2 



dy vr 2 7 ( x — z ) 2 (y _ z ) 2 

x [JV(x - y, b, 0) - N(x - z, h, 0) - N(y - z, b, 0) + iV(x - z, b, 0)N(y - z, b, 0)] 2 . 

III. INCLUDING QUANTUM EVOLUTION 

When the collision energy becomes high enough, multiple gluon emission becomes possible. 

Parametrically, each gluon emission brings in a factor a s ln(l/x). Accordingly, quantum evolution 

i_ 

takes place when x < e a « . Let the incident onium be characterized by the two-vector r. In 
course of evolution dipoles of different sizes are produced until eventually a dipole of size r' emits 
a gluon at rapidity y with transverse momentum k. In terms of the Regge theory, evolution in the 



rapidity interval between the original onium and the emittec 
single Pomeron, in agreement with the AGK cutting rules 27|, 



gluon corresponds to exchange of a 



65]. Afterwards, i.e. in the rapidity 



interval between the emitted gluon and the target nucleus, evolution is non-linear and corresponds 
to exchange of diffractively cut fan diagram, see Fig. [2] and Fig. [3l The general method for including 
the non-linear low-x evolution into inclusive processes in the dipole model framework was derived 
in [2^] and is applied later in this section. 

We would like to separately consider the following two cases: (i) rapidity gap Yb equals the 
produced gluon rapidity y, see Fig. [2j and (ii) a more general case Yq < y, see Fig. [3l In later 
sections we will focus our attention on the former case. The reason is that experimentally, diffractive 
production is usually measured per unit of invariant mass of the diffractively produced system 
(rather than per dk 2 ). The invariant mass is given by M 2 = k 2 /x where k and x = e~^ ¥ ^ y ^ refer 



5 



to the slowest particle in the gluon cascade originating from proton. That being the case, it is 
sufficient to consider production of a gluon adjacent to the rapidity gap as depicted in Fig. [2J This 
case corresponds to the rapidity gap Yq being equal the rapidity of the produced gluon Yq = y. 



A. Gluon production with y = Yq 




y 



FIG. 2: Fan diagram describing the diffractive gluon production in onium (P) - heavy nucleus (A) collisions 
with the rapidity gap being equal to the rapidity of the produced gluon. 



The relevant fan diagram is displayed in Fig. [2j We include the evolution effects using the 
method derived in [27]. We obtain the following generalization of (J6j) : 

da(k,y) _ a s C F I f , 2 , , 2 _ f , 2 , _ , , _ .„ , ._, , _ , |2 



d 2 k dy 



"f 2 F (2 ^ )2 J d 2 bd 2 B J d 2 r'nx{v,v\B-h,Y- y )\l{r'^ 



(7) 



where 

I(x'-y',k,y)= / d 2 z ie -^ 



zi — x 



Zl 



,r\2 



rl\1 



x [TV (x ; - y', b, y) - iV(x' - Zl) b, y) - N(y' - z u b, y) + 7V(x' - z 1; b, y)N(y' -z u b, y)] (8) 

Here ni(x — y, x' — y', B — b, Y — y) has the meaning of the number of dipoles of size x' — y' 
at rapidity Y — y and impact parameter b generated by evolution from the original dip ole x — y 
having rapidity Y and impact parameter B 63]. It satisfies the BFKL equation 
9m(x-y,x'-y',b,y) _ a s N c [ ^ (x - y) 2 



66 



67] 



dy 2tt 2 J (x — z) 2 (y — z) 2 

[ni(x - z,x' - y',b,y) +m(y - z,x -y',b,y) - m(x - y,x' - y',b,y)] , (9) 
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with the initial condition 

m(r,r',b ) 0) = tf(r-r')$(b), (10) 

where we denoted r = x — y and r' = x' — y'. 

The forward elastic dipole-nucleus scattering amplitude satisfies the nonlinear BK equation 



68, 



69] 



dN(x -y,b,y) _ a s N c [ (x - y) 2 

dy 2vr 2 

+iV(y - z, b, y) - iV(x - y, b, y) - iV(x - z, b, y)N(y - z, b, y)] , (11) 



with the initial condition given by ©. In writing both equations ([9]) and (jlip we assumed that 
the absolute value of impact parameter b is much larger than the typical dipole size. This is a 
justified approximation for a scattering off a heavy nucleus. 

In order to keep expressions as compact as possible, it is convenient to assume that the nuclear 
profile is cylindrical. This simple model allows correct identification of the atomic number (i.e. 
centrality) dependence of the cross sections. An explicit impact parameter dependence, which is 
required for numerical analysis, can be easily restored in the final expressions. Since we are not 
concerned here with the details of the impact parameter dependence, we integrate §§§ over b. The 
quantity 

n p {r, r',y) = J d 2 b n p (r, r', b, y) (12) 
in turn satisfies the BFKL equation with the initial condition 

n p (r,r',0) =8(r-r'). (13) 

In terms of n p (r,r',y), (J7|) reads 

da(k,y) a s CF 1 



■S A J d 2 r'n p (r,r',Y-y)\I(r',k,y)\ 2 , (14) 



d 2 kdy tt 2 (2tt) 2 
where Sa is the cross sectional area of the interaction region. 

The total cross section for diffractive gluon production is convenient to write in the following 
form 

da(y)_ = <^sCf Sa f ^t' n p (r,r\Y - y) [ d 2 w ^— ■= — ■= 

dy tt z J J (w — r y w z 

x [N(r', b, y) - N(w - r', b, y) - N(w, b, y) + N(w - r', b, y)N(w, b, y)] 2 . (15) 

where we introduced a new variable w = z — y' such that w is the size of one of the daughter 
dipoles formed by emission of a gluon at point z by a parent dipole r' = x' — y'. 
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B. Diffractive production with Yq < y 




So far we have been concentrating on a case in which the rapidity of the produced gluon y 
coincides with the rapidity gap Yq in a diffractive event. In this case the diffractive scattering 
amplitude Nd(t, h,y,Yo) coincides with the square of the forward elastic scattering amplitude 
N(r,b,y). In principle, a question may arise about the diffractive production of a gluon with 
y > Yq. Such process is shown in Fig. [3j For the processes in which the transverse coordinate of 
the gluon in the amplitude zi is approximately the same as its coordinate in the c.c amplitude Z2, 
the corresponding cross section is still given by (|14j) 

da{k,y) u s Cf 1 



;S A J d 2 r' n p (r,r> ,Y - y)\I(r> ,k,y;Y )\ 2 , (16) 



d?kdy tt 2 (2tt) 2 
where now in place of ([8]) we write 

I(x -y,k, y; y ) = y <iz ie ^^-—^-——-j 
- JV| (x' - zi, b, y; y ) - iV| (y' - zi, b, y; Y ) + iV| (x' - Zl , b, y; F )iV| (y' - Zl , b, y; F ) 



^( x '-y',b,y;y 



(17) 



The amplitude Nd(t, b, y; Yo) equals to the cross section of single diffractive dissociation of a 
dipole of transverse size r, rapidity y and impact parameter b on a target nucleus. It satisfies the 



Kovchegov-Levin evolution equation [7 1 
dN D (x-y,h,y;Y ) 



dy 

2ql s Cf f , 2 
= — — ' " z 

a sCF I ,2 



(x - y) : 



27r<J(y-z)ln(|x-y|A) 



7T^ 



(x - z) 2 (y - z) 2 

.2 

[N D (x - z, b, y; Y )N D {y - z, b, y; Y 



N D (x- z,b,y;Y ) 



(x - z) 2 (y - z) 2 



4iV D (x - z, b, y; Y )N(y - z, b, y) + 2iV(x - z, b, y)iV(y - z, b, ; 



(18) 



with the initial condition 



N D (r,b,y = Y ;Y ) = N 2 (r,b,Y ) . 



(19) 



Diffractive gluon production of the kind shown in Fig. [3] requires a dedicated study and will 
certainly lead to a number of interesting observations. We are going to perform such analysis in 
future publications. In this paper we concentrate on the case y = Yq. 



IV. DIPOLE EVOLUTION IN ONIUM AND NUCLEUS 
A. Dipole evolution in onium 

Dipole evolution in onium is encoded in the function ni(r, r', b, y) and is determined by solving 
equation © with the initial condition (jlOp . In the case of a cylindrical profile we use function 
n p (r,r',y) instead. The general solution of the BFKL equation for n p (r,r',y) reads 

/oo 
due 2 ^ x ^ y (r/r') 1+2iv CP, (20) 
-oo 

where a s = a s N c /TT and the leading BFKL eigenvalue 

X (u) = V(l) - \^{\ ~ ™) ~ \^{\ + w) , (21) 
with ip(v) being the digamma function 

= (22) 
The Mellin image Cf, can be found using the formula 

S(r-r') = ^ J_Jv{r/r')^. (23) 

The result is 

n P^ r '' V) = ^Fi J dv e2&sX{V)V (r/r') 1+2 ™ . (24) 
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Integral over v can be done analytically in two important limits. In the leading logarithmic 
approximation (LLA) we expand the function xiy) near the minimum at v = as 

x(%LA~21n2-7C(3)i/ 2 , (25) 

where Q{z) is the Riemann zeta function. Substituting (I25p into (124ft and integrating around the 
saddle point 

_ iln(r/r') 

" 14C(3)a sy ' (26) 

we derive 



n P (r,r',y) LLA « /IZZ3 e (^-i)y e wHi, a s y » ln 2 (r/r') , (27) 

where ap — 1 = 4a<j In 2. 

Alternatively, we can expand xiy) near one of its two symmetric poles at 2iv = ±1. This 
corresponds to the double logarithmic approximation. The choice of a particular pole depends on 
the relation between r and r' . Expanding near 2iv = 1 we attain 

X(-W-T^-- (^) 



Plugging this into (|24|) we have 

Mr, r', u)dla *^r 2 fjv e^ + ^ ^ . (29) 
The saddle point of the expression in the exponent is 



^('-Vh^J- (30) 

which is valid only if r < r'. Expanding the argument of the exponential near the saddle point u* 
to the second order and integrating gives the double-logarithmic approximation 

n P (r,r',y)nLA « J^^lf^ e V=W'/r) , r</ , H r'/r)»a s y. (31) 

To derive an analogous expression at r > r' we expand xi v ) near the symmetric pole 

(32) 

In analogy to ([30]) and (j3"T|) we derive the saddle point 



and the dipole density 

^r'.W^^e 2 ^, r>/, ln(r/r') » a s y . 



(34) 
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B. Dipole evolution in a heavy nucleus 

1. Leading twist approximation 

Consider the forward elastic dipole-nucleus scattering amplitude 2V(r, b, Y) satisfying the non- 
linear evolution equation (jlip . If the dipole size is much smaller than the saturation scale Q s , then 
the quantum evolution of the amplitude is governed by the BFKL equation. Therefore in this case, 
the general solution is 

/oo 
dve m ^y(rQ s0 ) 1+2 ™C^. (35) 
-oo 

The Mellin image of the amplitude N(r,h, 0)lt is calculated as follows 



c a = Qso r dr{rQso) -2-^ NirMO)LT 

n Jo 



QsO r, (n ,-2-2, 1 2 n 2 1 l + (l-2w)ln^f 

/ ar rCJ,n -r (J„ n = ; rn — — 

vr J y v u; 8 Vs0 8vr (1 - 2iv) 2 



In the last line of (j36j) we used the fact that Q s q logarithmically depends on r, see ((3|),(j4]). Analo- 
gously to the derivation of (|27p we obtain in the leading logarithmic approximation 



(37) 



where the saddle point is 



In the double logarithmic approximation ([28]) the saddle point for the case r < 1/Q s q is 



»V s o 



Repeating the by now familiar procedure we write 

ln 1 / 4 ^ 



N(r,b,y) DLA - i67T {2 _ sy)3/4 rQ^yi + ^J In ^ jey 



r < 1/Q s0 , In — — > a s y . (40) 
rQ. s o 



Next, we consider the case r > l/Q s 
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2. Deep saturation region 



71 



Q,EJ 



Solution to the BK equation (jlip deeply in the saturation regime was found in 
With the logarithmic accuracy the dominant dipole splitting corresponds to the configuration in 
which the size of one of the daughter dipoles (~ 1/Q S ) is much smaller than the other (see Sec. IVip . 
Denote again r = x — y and w = z — y. In the saturation region we have either w <r ~ |w — r| 
or the symmetric configuration |w — r| <C r ~ uu. Both give equal contribution to the integral over 
w. Restricting ourself to the case w <C r and doubling the integral we write the BK equation as 
follows: 

dN(r,b,y) a s Cp f r2 dw 2 r „ T , . „ T , . , T/ .„ , 

*^-2 — ^ [N ( w > b, y) — N(w, b, y)N(r, b, y)] . 41 

dy vr J 1/Q 2 w 2 

Now, for the reason that in the saturation region, the amplitude N(r, b, y) is close to unity we 
render (14T1) as 



d{l- N(r,h,y)} a s C F [ dw 2 2a s C F 2 2 

- 2/ _ _{i_JV( r b,y)} = ln(r QJ{1 - N(r,b,y)}. 

(42) 

The saturation scale Q x {y) c an be found by equating the argument of the exponent in (f4H to a 



<9y 



constant which yields 



71 



7| 



Qs(y) ~ Qsoe 2 "^ • (43) 



r energy limit of the 



Introducing a new scaling variable r = ln(r 2 Q^) we solve (j42]) and find the hig 
forward scattering amplitude (in the fixed coupling approximation). It reads 

N(t, b,y) = l-S e" r2 / 8 = 1 - 5 e"i ln2 ( r2 «?) , (44) 

where we approximated ~ -/V c /2 in the large iV c limit. S'o is the integration constant. It 
determines the value of the amplitude at the critical line r(y) = 1/Q s (y). 

V. DIFFRACTIVE CROSS SECTION IN THE QUASI-CLASSICAL APPROXIMATION 

Careful inspection of JT]) reveals that the cross section vanishes when size of the onium r = x — y 
is much larger than the characteristic scale 1/Q s q. This is in a sharp contrast with the inclusive 



gluon production case [26| where the cross section stays finite at r — ► oo. To understand the reason 
for such different behavior, consider a sample diagrams contributing to each of the processes shown 
in Fig. [H In diagram (a), corresponding to the inclusive case, the propagator of the qqg system 
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y 



y 



p 



p 




X 




X 



(a) 



(b) 



FIG. 4: An example of diagrams contributing to (a) inclusive gluon production and (b) diffractive gluon 
production. 



contrary, in diagram (b), corresponding to the diffractive case, the propagator involves both x 
and y coordinates, see ([2]), no factorization of y dependence similar to the inclusive case happens. 
All other diagrams contributing to these two processes can be analyzed in the same way. Let us 
look at the color structure of the qqg system in the two cases. By Pomeranchuk theorem, the 
exchanged (Coulomb) gluons are in the color singlet state. Therefore, we notice that in diagram 
(a) the qqg system is in the color octet state and its propagator through the nucleus equals the 
propagator of a gluon dipole, whereas in diagram (b) the qqg system is always in the color singlet 
state corresponding to the quark dipole. This feature can be seen also in expressions for the 
propagator, one involving the gluon saturation scale Q 2 sQ (inclusive case), another involving the 
quark saturation scale jf-Q% ~ \Qlo (diffractive case). 

Now we would like to determine how does the diffractive cross section behave in the quasi- 
classical approximation. Consider the following integral appearing in the r.h.s. of (|7|) and (|15|) : 



Let us analyze its behavior in the quasi-classical case (y = 0) for small r < l/Q s o and large 
r > l/Qso onium sizes. 





iV(w - r, b, y) - N(w, b, y) + N(w - r, b, y)N(w, b, y)} 2 . 



(45) 
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A. Dilute regime r < l/Q 8 o 

In the case of small onium we divide the entire integral over w into three terms as follows: 

„2 



J(r,0) « / -4- [l-iV(r,b,0)] 2 iV 2 (w,b,0) 

+ r 2 / ^ [^fo b > °) " 2iV ( w > b > °) + ^ ( w > b ' °)] 

J r ^ 

+ r 2 / ^- [N(r, b, 0) - 27V(w, b, 0) + iV 2 (w, b, 0)] , (46) 
Ji/Qso w 

where (|46p holds in logarithmic approximation. We can estimate each term utilizing the fact that 
according to ([5]) 

{1 , r > I/Qsq , 

(47) 
\r 2 Q 2 sQ , r«l/Q s0 . 

Moreover, in the last two terms in (|47p N(r, b, 0) can be neglected since in most of the integration 
regions w S> r. Indeed, in the second term in the r.h.s. of (|46p it can be seen once we neglect 
iV 2 (w,b,0) and expand N(w, h, 0) at small w. In the third term in the r.h.s. of (|46p we have 
2A r (w,b,0) — iV 2 (w,b,0) ~ 1 while 2V(r,b, 0) <C 1. Accordingly, expanding the integrands using 
(I47p we find that the first term in the r.h.s. of (|46f) is of the order O(r s Q^ ), whereas the second 
and the third ones are of the order O(r 2 Q 2 ). Therefore, the last two terms in (I46p dominate in 
the regime rQ s o C 1. It is customary to denote 

N G (r, b, y) = 2N(r, b, y) - N 2 (r, b, y) , (48) 

which has the meaning of the gluon dipole forward elastic scattering amplitude. In terms of this 
quantity the function J(r, 0) reads 

f°° rim 2 

J(r, 0) « r 2 / ^ iV 2 (w, b, 0) , r « 1/Q a0 , (49) 
Jo w 

where the lower limit of integration (r) has been set to zero with logarithmic accuracy (note that 
the second integral in the r.h.s. of (|46p is dominated by dipoles of size w ~ l/Q s o S> r). 

It is useful to notice, that (|49p holds also in the case the low-x evolution is taken into account. In 
the quasi-classical approximation the integral (|49p can be done if we treat Q s q as a w-independent 
constant neglecting its logarithmic variation. In that case, substituting ([5]) we derive 

J(r,0)«^ln2r 2 Q 2 , r « l/Q s0 , (50) 

The cross section is obtained using © and ([431) . We have 

da o> s Cf a j i n x a s Cpln2 2 2 

— = 5 , j4 J(r,0) = Sat Q s0 , r«l/Q s0 . (51) 
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B. Dense regime r > 1/Q 



As in the previous case we divide the integral into three parts 

•l/Q s o dw 2 



J(r,0) « / [l-iV(r,b,0)] 2 iV 2 (w,b,0) 

Jo ^ 

f ^ [l-7V(r,b,0)] 2 iV 2 (w,b,0) 
Jl/Qso w 

f°° dw 2 9 
+ r 2 / — [JV(r, b, 0) - 2N(w, b, 0) + iV 2 (w, b , 0)] . (52) 

J r ^ 

Utilizing (H7|) we simplify (152]) in the logarithmic approximation as follows 

/c, o) - w ,, b, o, - 1]» ( r-° g + r £ + h r d 4) ■ m 

\/o ^ 2 64 7 1/QsQ 7 r u; 4 J 

« [iV(r,b,0)-l] 2 ln(r 2 Q 2 ), r » 1/Q s0 . (54) 



Going from (|53p to (|54p we kept only the second term in the brackets in (|53p as it is logarithmically 
enhanced. Formula (|54p is valid in the case of low-x evolution as well. Substituting ([5]) into (|54p 
yields 

J(r) = ln(r 2 Q 2 ) e^^o , r » i/Q s0 . (55) 
Finally, the cross section follows from ([7|. (|45jl and (|55p as 

^ = ^^ln(r 2 g 2 )e-^o, r » 1/QsQ . (56 ) 

The striking feature of this formula is strong exponential suppression of diffractive gluon production 
for large onium. We will see in the next section that this result completely changes when the 
quantum evolution in the onium becomes an important effect. 



VI. DIFFRACTIVE CROSS SECTION INCLUDING LOW-i EVOLUTION 



Using (|15p and (145)) we write 



do _ OsCf_ ^ f d 2 r , ^ r / y — y)J (r', y) . (57) 
dy n J 
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A. Dilute regime r < 1/Q s (y) 

As in the quasi-classical case, first we are going to find the kinematic region which gives the 
largest (logarithmic) contribution to the integral. We have 

da o s Cf 



dy 7r 



S a 2tt 



dr r n p (r, r , Y — y) J(r' , y) 



/ dr'r'n p (r, r' , Y — y) J(r', y) + / dr'r'n p (r, r' , Y — y) J(r', y) 
Jr Ji/Q, 



(58) 



As has been noted in the previous sections, equations (|49p and (|54p hold also in the evolution 
case, provided the y-dependence is explicitly indicated in the arguments of J(r,y) and iV(r, h,y). 
Generalization of (|50p reads 

J(r',y)^C r' 2 Q 2 (y), r' < 1/Q s (y) , (59) 

where Co is a constant which depends on a particular functional form of Nc(r,h,y) and can be 
found numerically from (|11|) , Using (|44|) in (|54|) gives another limit of function J(r',y): 

J(r',y) « <S 2 e- ln V<W ln(r' 2 Q 2 s ) , r' » 1/Q s (y) . (60) 

Accordingly, using (|3ip or (|34|) depending on the relation between r and r', i. e. n p ~ r 2 / r /4 if 
r < r' or n p ~ 1/r' 2 if r > r', as well as (|59p and (|60p we estimate that the second integral in (|58p 
is enhanced by In -k- with respect to the first one, whereas the third integral is vanishingly small. 
Thus, 

^ ^^L SA 2n [ 1/QS dr'r'n p (r,r',Y-y)J(r',y) 
dy tt J r 

_ 2C a s C^ A , 2 rVO. (2« s (y-y))V 4 2 ^ s( y- s)b r! 

47r 3/2 ^ I r , ^3/4^ • [pl) 

Changing to a new integration variable rj defined as rj 2 = In — the integral in (|6ip can be taken 
explicitly in terms of the imaginary error function. In the double-logarithmic approximation the 
result reads 

T ~ C °Hf A <■>- IV - W 4 ^ m - (Y -"""^ ■ r «!/<«•<,). (62) 

dy An^ 2 (2a s {Y - y) ln^) 1 / 4 

Both the quasi-classical result (|5ip and its quantum counterpart (|62p show that the cross section 
is proportional to r 2 as required by the color transparency. 
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B. Dense regime r > 1/Q S 

Analogously to (|58|) we get 



da a s C F 

■ DA ^ 



dr'rUpiv, r', Y — y) J(r' , y) 



dy ir 

rr poo 

+ / dr'r'n p (r, r', Y — y) J(r', y) + / dr'r'n p (r,r' ,Y — y) J(r' , y) 
■/i/Q s -A- 



(63) 



The logarithmically enhanced contribution arises from the second integral which - upon substitu- 
tion of (l60l) and (l34l) - becomes 



ckr _ QsCf^a „ 2 ( (y _ vv V 4 T <fr' ln(r /2 Qg) i n(r 'Q s) 2^(1^)1^ ffi ^ 



r 



Note that in the relevant kinematic region 1/Q S « r' C r we can approximate In p- = ln(rQ s ) + 
In ~ ln(rQ s ). The integral over r' then becomes trivial yielding the final result 

T = ^0^0 {m : [ L;f' A eV^(^)»n(r Q .) , r >l/0. (y ). (65) 

27T 3 / 2 ln 3/4 (rQ s ) 

We observe that the cross section given by ()65[) is an increasing function of the rapidity interval 
Y — y between the onium and the nucleus. Together with the quasi-classical expression (|56p it 
implies that the total cross section for the diffractive gluon production in a scattering of a large 
onium off the heavy nucleus is non-vanishing only if the low- a; evolution in onuim is an important 
effect. This can be seen directly in Fig. [2j the gluon multiplicity arises from the cut Pomeron 
attached to the onium. This observation has important phenomenological consequences as we 
discuss in the next section. 

VII. DISCUSSION AND SUMMARY 

In this paper we discussed the coherent diffractive gluon production in high energy onium- 
nucleus collisions. The gluon multiplicity in the case of onium of small size r < 1/Q S is given by 
(|5ip and (|62p and can be summarized as follows 

oc r 2 Q 2 s {y) xG (exp(y - Y), Q 2 S ) , r « Q s , (66) 

where xG(x,Q 2 ) is a gluon distribution function at momentum scale Q 2 . Gluon multiplicity van- 
ishes in the limit r — > as is required by the color transparency. 

In the other limit of large onium, the gluon production cross section vanishes in the quasi- 
classical approximation as implied by ([56|) . At a s {Y — y) > 1 the evolution effects in onium 
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play increasingly important role. It is the cut Pomeron, connecting the onium and the dipole (r') 
emitting the triggered gluon, which contributes to the fast increase in gluon multiplicity as interval 
Y — y increases. One way to see it is to recall that during the linear evolution dipoles of various 
sizes are produced from the parent onium of size r. We explained in (|63p and (|64p that the main 
contribution to the multiplicity stems from the dipoles of size r' ~ 1 /Q s no matter how big is the 
initial dipole r. The resulting expression (|65[) has the following behavior 

dN D (y) 



OC X 



G(exp(y -Y),Qf) , r > Q s , (Y - y) > l/a s . (67) 



dy 

Dependence of the diffractive gluon multiplicity on the onium size is summarized in the Fig. [5) 



dN 1 
dy 




FIG. 5: Relative multiplicity of diffractive gluons as a function of the onium size r in the quasi-classical 
approximation (labeled y = 0) and at very low x (y = 2). The rapidity values correspond to those at RHIC 
as explained in the text. We chose Q s o = 1 GeV and Q s = 1.35 GeV. iVoo is a normalization constant. 



To the extent that the large onium can serve as a model for proton, (|56p and (|65p describe the 
diffractive gluon production in proton-nucleus collisions. As such it has a direct implications to 
the RHIC and LHC phenomenology. Phenomenological studied show that the gluon saturation at 
RHIC starts to impact the gluon and valence quark spectra at rapidities 77 ~ 1 (and larger). In our 
notations it corresponds to the rapidity interval y ~ 6 between the gluon and the heavy nucleus 
and Y — y ~ 4 between the proton and the gluon. This corresponds to x p ~ e -4 ~ 0.02 which 
is perhaps insufficient to have a sizable low-x effect in proton implying a very low multiplicity of 
diffractive gluon production. On the other hand, exploring the backward rapidity region r/ < 
will not allow to probe the gluon saturation in the nucleus. Therefore, if the typical inter-quark 
distance in proton is larger than ~ 1/Q S — 0.2 fm, then we do not expect a significant multiplicity 
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of gluons in coherent diffraction of a proton on nucleus at RHIC, see Fig. [5j 1 

The situation radically changes at LHC where an additional rapidity window Ar/ ~ 6 opens up. 
From the point of view of gluon saturation, the mid-rapidity in pA at LHC is expected to be similar 
to the rapidity r\ = 3 at RHIC [^Q]- At the same time, at the LHC midrapidity, x p ~ e~ 7 = 0.001 
which is certainly sufficient for the low-x evolution to take place in proton. Therefore, we expect 
that measurements of the diffractive gluon production in pA collisions at LHC will be a sensitive 
probe of the low-x dynamics. At EIC the typical dipole size r is determined by the photon virtuality 
Q as r ~ 1/Q which makes it possible the detailed study 
gluon production in different kinematical regions (see e.g. [75l . 

An interesting extension of our work is a case of diffractive production with small rapidity 
gap Yq < y, which will be relevant at LHC. Of special interest is dependence of the differential 
cross section for diffractive gluon production on transverse momentum of produced gluon. We are 
addressing this and other issues in the forthcoming publication. 



uon saturation using the diffractive 

,Q). 
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